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Abstract. We present a simple technique to compute moments of derivatives of unitary characteristic polynomials. 
The first part of the technique relies on an idea of Bump and Gamburd: it uses orthonormality of Schur functions over 
unitary groups to compute matrix averages of characteristic polynomials. In order to consider derivatives of those 
polynomials, we here need the added strength of the Generalized Binomial Theorem of Okounkov and Olshanski. 
This result is very natural as it provides coefficients for the Taylor expansions of Schur functions, in terms of shifted 
Schur functions. The answer is finally given as a sum over partitions of functions of the contents. One can also obtain 
alternative expressions involving hypergeometric functions of matrix arguments. 

Resume. Nous introduisons une nouvelle technique, en deux parties, pour calculer les moments de derivees de 
polynomes caracteristiques. La premiere etape repose sur une idee de Bump et Gamburd et utilise l'orthonormalite 
des fonctions de Schur sur les groupes unitaires pour calculer des moyennes de polynomes caracteristiques de matri- 
ces aleatoires. La deuxieme etape, qui est necessaire pour passer aux derivees, utilise une generalisation du theoreme 
binomial due a Okounkov et Olshanski. Ce theoreme livre les coefficients des series de Taylor pour les fonctions de 
Schur sous la forme de "shifted Schur functions". La reponse finale est donnee sous forme de somme sur les parti- 
tions de fonctions des contenus. Nous obtenons aussi d'autres expressions en terme de fonctions hypergeometriques 
d' argument matriciel. 

Keywords: random matrix theory, hook-content formula, moment of characteristic polynomials, shifted Schur func- 
tion, generalized Pochhammer symbol, hypergeometric function of matrix argument 



1 Introduction 

We take for the characteristic polynomial of a TV x N unitary matrix U 

N 

Zu(0) := n(l-e^- e )), (1) 

where the OjS are the eigenangles of U and set 

Vu (0) := e iN ^/ 2 e~ l ^ e ^ 2 Zu(0). (2) 
It is easily checked that for real 9, Vtj (6) is real and that Vjj (0) equals | Zjj (0)\. 
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For k and r integers, with < r < 2k, we will investigate the averages (with respect to Haar measure) 



(M) N (2k,r) 
\V\ N (2k,r) 



Zu(0) 



\Vu(0) 



2 k 



2 k 



\Zu(0)J / U(JV) 



Vu(p) 



(3) 
(4) 



XJ(N) 



This is notation we already used in [Deh08] (where a |7W was also present but is not needed here), and 
is only notation in the LHS: the (Ai) and |V| are thus meant each as one symbol and are supposed to 
mnemotechnically remind the reader of what is in the RHS. We immediately state the following easy 
lemma. 



Lemma 1 For k > h non-negative integers, we have the relation 

2h 



\V\ N (2k,2h) 



E 

i=0 



2h 



(M) N (2k,i) 



2ft- 



(5) 



Proof: This is available (in the same notation) in [ Deh08 ] and a consequence of Equation (|2j), which leads 
to the polynomial relations between and their norms. These relations give 



\V\ N (2k,2h) = 



I'M' ^V' ; l( J 



3=0 s ' 1=0 

which is easily deduced from [ peh()8| |. 

We are actually more concerned with the renormalizations 

(M) N (2k,r) 

]\fk 2 +r 

V\ N (2k,r) 



(M) N (2k,j + l), 



(M)(2k,r) = lim 



\V\(2k,r) 



lim 

7V->oc 



]\fk 2 +r 



(6) 

□ 

(7) 
(8) 



Theorem g will show that these normalizations are appropriate. 

The random matrix theory problem of evaluating (M) (2k, r) and | V| (2k, r) has applications in num- 
ber theory (see [ |Deh08 1 for a more detailed exposition of these ideas). Indeed, these values are related to 
the factor g(k,h) in the formula 



lim 

T^oo T (log 



T_\k 2 +2h 



2k-2h 



2h 



At 



a(k)g(k,h), 



(9) 



where a(k) is a (known) factor defined as a product over primes. This, along with a discrete moment 
version due to Hughes, is the principal underlying motivation for the all the random matrix theory analysis 
that occurs in jHKO00| , ftug01| , ^ug05j |Mez03[ pRS06| , |FW0^ |. 



Moments of derivatives of characteristic polynomials 



3 



Another application is tied to the work of Hall frlal02a Hal02b, Hal04, |Hal08 ], where results on the 
objects studied here can be used to hint towards optimizations of rigorous arguments in number theory, and 
serve as (conjectural) inputs on theorems there. The number theory statements concern average spacings 
between zeroes of the Riemann zeta function. The works of Steuding [ |Ste05 ] and Saker [3ak09] follow 
similar approaches. 



In [Deh08|, the author investigated for fixed r ratios of those quantities and established they were a 



rational function. We reprove this result here, but with a much simpler method leading to a much simpler 
result. In particular, rationality of the RHS is transparent from the following statement (definitions are 
given in Section ||), since the RHS sums are finite sums over partitions of r: 



Theorem 2 For < r < 2k, with r, k £ 



(M) N (2k,r)_. r 



{M) N (2k,0) 
(M)(2k,r) :7 



v H (JVt /*)((-*) tM) 

fcK (-2*0 t/i ' 

r\ k t M 



(M)(2k,Q) 

while the denominators on the left are known ( [ BG06 /, flfiSOO J ): 

G(N + 2k + 1)G{N + 1) G(k + l) 2 



(10) 
(11) 



(M) N (2k,0) 



{1} 



2A- 



G{N + k + iy 



G(2k + 1) 



(12) 



and 



(M)(2k,0) 



where G(-) is the Barnes G-f unction. 



lim 

N->oc 



(N k ) 



({I} 2 " 



N k2 



G{k + 1) 2 
G{2k + lY 



(13) 



We now briefly discuss the technique used to obtain this theorem. The first idea will be similar to an 



idea of Bump and Gamburd [BG06| of using orthonormality of Schur functions to efficiently compute 
The new idea here is to combine this with the Generalized Binomial Theorem (E3|) 



matrix averages. The new idea h ere is to combine this with the Generalized Binomial Theorem (f23p of 
Okounkov and Olshanski [ OQ97 1 in order to obtain information about the moments of derivatives instead 
of moments of polynomials directly. 

This paper is structured as follows. We give in Section || the basic definitions needed. In Section ^, 
we explain the Generalized Binomial Theorem. We prove Theorem || in Section [I| We use this result 
to deduce in Section |] further properties of the rational functions obtained. We present in Section ^ 
an alternative interpretation of these results in terms of hypergeometric functions of a matrix argument. 
Finally, we announce briefly in Section ^[further results. 



2 Definitions 

Since Theorem || presents its result as a sum over partitions, we first need to define some classical objects 
associated to them. We follow conventions of [Sta99| throughout. 
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Fig. 1: The Ferrers diagram of partition (9, 6, 2, 1). 
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Fig. 2: The contents and hook lengths of the partition (9, 6, 2, 1). Its hook number is thus 4180377600. 

Partitions are weakly decreasing sequences Ai > A2 > • • ■ \(\) of positive integers, its parts. The 
integer I = 2(A) is called the length of the partition A. We call the sum ^\ Ai of its parts the size | A| of the 
partition A. We sometimes say that A partitions |A|, which is written A h |A|. If the partition has k parts 
of equal size N, we simplify notation to (N k ). 

We always prefer to think of partitions graphically. To each partition we associate a Ferrers diagram, 
i.e. the Young diagram of the partition presented in the English convention. We only give one example 
(Figure [j]) as it should be clear from it how the diagram is constructed: the parts A* indicate how many 
standard boxes to consider on each row. 

There exists an involution acting on partitions, which we denote by A 4 . Its action on diagrams amounts 
to a reflection along the main diagonal. 

Partitions can be indexed in many different ways. Indeed, we have already seen that finite sequences 
of (decreasing) part sizes can be used as an index set. Shifted part lengths are essential for the work of 
Okounkov and Olshanski underlying Section^, but we do not need to define that system of coordinates. 

Define the content of a box □ located at position (i, j) in a partition A as c(D) = j — i (see Figure [|). 
We use this to define the symbol 

ktn := H(k + c(n)). (14) 

This definition leads to k t (n) = k(k + 1) • ■ • (k + n - 1) and k t (1") = k(k - 1) • ■ • (k - n + 1). 
We sometimes abbreviate the first k ^ n and the second k 4- n. This is clearly a generalization of the 
Pochhammer symbol. Indeed, we even adopt the convention that N f (—k) = l/((N + 1) t k), which 
guarantees 

{N + a- 1) I (a + b) = (N t a) ■ ((N - 1) 4, b). (15) 
To get back to the generalization to partitions, we have an immediate relation under conjugation: 

fctM* = (16) 
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Given a box □ £ A, define its hook (set) 



{(i,j')&\:j'>j}U{(i',j)e\:i'>i} 



(17) 



Remark that the box □ itself is in its hook. Define the hook length \huj~, | as the cardinality of the hook 
(see Figure 0) and call their product the hook number of a partition A: 



hx := 



IW 

□eA 



(18) 



Hook numbers are of importance thanks to the hook length formula of Frame, Robinson and Thrall 



[FRT54]. This counts the number fx of standard tableaux of shape A, which is also the dimension 
dimx A = X A (1) of the character associated to the partition A for the symmetric group Si\\ (see [SagOl 



Sta99|) 



h 



|A|! 



(19) 



The last classical combinatorial object we need is the Schur functions . To each partition A we associate 
an element Sa of degree |A| of the ring Aq [X] of polynomials symmetric in a countable set of variables 
X = {xi}. We refer the reader to [ Bum04 | for definitions, and only state a few properties. 

Let r G N, and take a partition A of size r. One can define a map from \n from V(N) to C in the 
following way: 



Xn(9) ~sx(9) :=SA(e"\< 



,e if>N , 0,0,0,. 



)• 



(20) 



with the e i0 J the eigenvalues of g e V(N). When /(A) > TAT, XjvCflO = °' but once N - 'W' ^ Xw(ff) 
become (different) irreducible characters of U(iV). We thus have the formula 

(xat>Xw)u(ao = ( s ^('): s /*('))u(JV) = 



1 if A = and N > |A|, 
otherwise. 



(21) 



The last formula we need concerns the evaluation of Schur polynomials, at repeated values of the argu- 
ments. Denote by {a} R the multiset consisting of the union of R copies of a and countably many copies 
of 0. The hook-content formula (see [Sta99, Bum04]), a consequence of the Weyl Dimension Formula, 
states then that 



*a {iy 



fct A 



(22) 



3 Shifted Schur Functions and Generalized Binomial Theorem 



The Generalized Binomial Theorem as formulated in [ |Q097i Theorem 5.1] can be interpreted as a Taylor 
expansion of the character x„ = Sa(") of U(n) around the identity Id nx „. It says explicitly that 



S A (1 +£!,•■• , 1+X n ) 

*a({1}") 



E 

l(fi)<.n 



S ^(^li'"' i A„)s /i (xi, • • ■ ,£„) 

nf/i 



(23) 
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where the s* are shifted Schur functions. Those were introduced by Okounkov and Olshanski in [ 3097] 
and need not be defined here as we only need their values on a very limited set of arguments. This is 
given by the following lemma, which is very elegant and seems to be new. Note the (almost-)symmetry 
between k and N. 



Lemma 3 



*;({N} k ) = vv({i}") •*/*({!} 



|H l (HV)tM)(fetM) 



(24) 
(25) 



Proof: We will need two equations from [0O97|. Equation (1 1.28) tells us that 



s*(xi,->- ,x„) = det [h*._ i+j (xi +j - 1,- • • ,x n +j - l)]f^ 



(26) 



for large enough R (it is then stable in R). Equation (1 1.22) from [ OQ97 1 deals precisely with those h*: 



K ({N}> 

Combining these equations, we obtain 



(Nlr)-t) r {1} 



dot 



((N + j-l)lQn-i+j))-fi 



tH-i+3 



{1}' 



dot [(N t j)((N - 1) ; (jm - i)) • Ki-i+i (i 1 }") 



R,R 
R,R 



J((7V - 1) I ( Mi -i)) JJ(/y ti) det f,^_^ {1} 



\j =1 



H((N + i-l)if,i)\ V ({1}' 



]jN-c(D)\ 5fi ({!}' 



^ V ({l} JV )5 /1 ({l} fc 



(27) 

(28) 
(29) 

(30) 

(31) 

(32) 
(33) 



The fourth line follows from Equation ( |l5|), t he fifth from reorganizing a product over rows into a product 
over boxes, and the sixth from Equation (|22|). □ 



We are now ready to launch into the proof of Theorem I 
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4 Proof of Theorem g 



The method of proof will be very similar to the technique presented in [BG06]. In particular, both the 
Cauchy identity [Bum04] 

\l + x iyj = ^s X t(xi)s\{yj), (34) 

where xi and y 3 are finite sets of variables, and the asymptotic orthonormality of the Schur functions will 
again play a crucial role. However, the power of their technique is now supplemented by the General- 
ized Binomial Theorem, which will provide for a dramatic simplification of the arguments and results in 
[ |Deh08p . 

Proof of Theorem ||: We have 

N 

= n - e ~ i0j i 1 - eWt ) (36) 

= (-l) N detUZu(0) (37) 

and thus 



Zu(0) = (-lf^detf/ Z v (0) k (38) 
= (-l) kN s {kN) (U)Zu(0) k . (39) 

We use the Cauchy Identity from Equation ([34|) to obtain 

Zu( ai ) ■ ■ ■ Zu(a r ) = ^s At (t/)s A (-e- iai ,--- ,-e- m ") (40) 

A 

= ^(-l)l A ls At (f/)s A ( e - lai ,--- ,e-^). (41) 

A 

To the first order in small a, we have e~ la « 1 — ia, so 

Z'uW = ^(-l) |A W(t/) d 1 ---d r s x (l-ia 1 ,--- ,l-tOr)| Oi= ... =ap=0 , (42) 

A 

where 9^ := 9 ai . 

Putting everything together, we obtain 



|Zt/(0)| 2 



1fw)' = < - 1, '*™ ) ^- 

]T(-l) |A| s At (^) 9 1 ---5 r s x ((2fc-r) x {1} U {1 - ia u ■ ■ ■ , 1 - iar})| Ol= ... =Or=0 - ( 43 > 
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Just as in the original proof of Bump and Gamburd, orthogonality of the Schur polynomials kills all terms 
in the sum but one (where A* equals (k N )) under averaging over U(iV). Hence this simplifies to 

ft • • • d r s {Nk) {(2k -r) x {l}U{l-iai,.-. ,l-M)| Ol= ... =Or=0 - < 44 ) 
This is the perfect opportunity to apply the Generalized Binomial Theorem. We wish to set n = 2k, and 

-idi for 1 < i < r 
for r + 1< i < 2k { ' 



in Equation ( |23[ ) to get 



({A^} fe ) ft •■•fts jU (ai,-- - ,a r )| ... =a =0 

^— (46) 

The (additional) restriction on fj, is obtained because of the derivatives: since the Schur functions are 
evaluated at a* = 0, and s M is of total degree in the a^, we must have /i| = r for something to survive 
ft • • • ft. 

When /i h r, we haveQ 

ft — ft-^K--- ,Or)| Ol= ... =ar=0 = <W<1'>) = dim X5 ltl| = 7^ ( 4? ) 

since derivation and multiplication by power sums are adjoint. 

Combined with Lemma (||), Equation ( [To] ) and Equation (|22|), this gives Equation (|l0|), which then 
quickly implies (|l~ij). □ 

5 Properties of the Rational Functions 



In the RHS of formulas ( |10| ) and (|llj), we have a sum of rational multiples of ratios of polynomials in k 
(and N), hence rational functions of k. We will now explain some of the properties of these functions, 
which are easily deduced from Equation ( |Tl| ) and the limiting version of Lemma [l]. 

Proposition 4 For a fixed r, h 6 N, there exists sequences of even polynomials X r , Y r and X 2 h such that 
for all k > r, 2h, 

(M)(2k,r) ( iV X r {2k) 

(M)(2k,0) \ 2 J Y r (2k)' ( } 

|V|(2fc,2/t) X 2h (2k) 

|V|(2fc,0) Y 2h (2k)' } 

(j ' Observe that higher derivatives require knowledge of values of symmetric group characters at other group elements than the 
identity. 
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and such that X r and Y r are of the same degree, monic, and with integer coefficients. We also have 
deg X2h < deg Y^h and that Xik has integer coefficients. 

This Proposition leads or relies on a few easy facts. 

Remarks. 

• All polynomials are even due to Equation (JT^). 

• We emphasize that we are not quite taking the simplest possible form of X r (2k), X r (2k) and 
Y r (2k) here. For some (presuma bly sma ll) rs, some spurious cancellations will occur. However, 
an explicit expression, given in [ fc)eh08 ], can be obtained for the Y r such that all the previous 



statements are true. Let us just say that we are taking for Y r the common denominator not of every 
term in the sum over partitions A, but of every (l-or-2-terms-)subsum over orbits of the involution 
on partitions, after all the simplifications of the type 2 k+2i = \ l ^ at occur m (2fc)tA ■ 

With this convention, the zeroes of Y r are exactly at the odd integers between 1 — r and r — 1. 

Since the squares of dimension of characters of a finite group G sum to the order of G, or due 
to the existence of the Robinson-Schensted-Knuth correspondence between permutations of n and 
pairs of Young tableaux of the same shape partitioning n, we have starting from Equation ( |l^ ) the 
identity 

En! 
¥ = 1, (50) 

which defines the Plancherel measure on partitions of n. The polynomial X r is monic thanks to 
this last identity. 

The sum appearing in Equati on (pTj) is a special case of a problem studied by Jonathan Novak in 



Equation (9.21) of his thesis |Nov09| 



6 Hypergeometric Functions of Matrix Arguments 

For completeness, we now discuss an alternative way to approach the expression in (|ic|). Hypergeometric 
functions of a TV x N matrix argument M are generalizations of hypergeometric functions of a complex 
variable. Originally, this extension is defined on multisets of complex numbers, which can then be seen as 
eigenvalues of a matrix using the trick of Equation (^0|). There is extensive literature on those f unction s, 



most of it tied to multivariate statistica l analysis. For a recent, accessible presentation, consider [FW08] 
For M a TV X N matrix, we define JDGR94 |GR9l| ] a function of M as follows: 



F(n h ■ M) y IK=i°itA s A (M) 
p F q {a u b 3l M) (51) 

where S\(M) is the Schur polynomial evaluated at the eigenvalues of M. This generalizes the classical 
hypergeometric functions of a complex variable z to (the multiset of eigenvalues of) a square matrix 
variable M, via the following substitutions: 
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the sum over integers is replaced by a sum over partitions; 

the generalized Pochhammer symbol replaces the rising factorial; 

the extra factorial that is always introduced for classical hypergeometric functions (by convention 
then) is replaced by a hook number; 

powers of z are replaced by Schur functions of the eigenvalues of M. 



• They admit integral representations, closely related to Selberg integrals (see [Kan93 1). 

By forming an exponential generating series of Equation < JTp| ) and with the help of Equation 
are able to obtain the (confluent) hypergeometric function 

^ {M) N (2k,r) (izY ^ 1 (iVtM)((-fe)tM) jA| „ , , „, T , , 

We can then substitute for the RHS of this last equation many different expressions: the theory of hy- 
pergeometric functions of matrix arguments also involves integral expressions, differential equations and 
recurrence relations. However, since we are interested in asymptotics of these expressions for large N, for 
which little theory is developed, this is unfortunately of no real use at the moment. Note though that the 
hypergeometric function that appears is special, as it is onl y evalu ated at scalar matrices. In that special 
case, N x N determinantal formulas have been developed [ |GR85 1 (but not in the confluent case). 



7 Further Work 

The occurrence of the Plancherel measure is not a coincidence. In fact, much can be derived from this, 
and this will be the basis of further work: we will prove in a subsequent paper that the leading term of 
X2h has coefficient jj^n and is of degree 2h lower than X-xh- 
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